Preprint typeset in JHEP style - HYPER VERSION 



PI-COSMO-73 



Stochastic tunneling for strongly non-Gaussian 
inflationary theories 



Andrew J. Tolley'^ and Mark Wyman'' 

Perimeter Institute for Theoretical Physics 
31 Caroline St. N 
Waterloo, ON N2L 2Y5 
Canada 



atolley@perimeterinstitute.ca ^ mwyman@perimeterinstitute.ca 



Abstract: We reconsider the dynamics of stochastic or thermal tunnehng in theories hke 
Dirac-Born-Infeld inflation that have non-minimal kinetic terms and, as a result, strongly 
non-Gaussian perturbations. We first describe a local description of the tunneling process 
which gives results consistent with the standard Hawking-Moss tunneling. This result is 
under perturbative control as long as the fluctuation determinant is well approximated 
by a one-loop integral. We then move to a global description, using the methodology of 
stochastic inflation and the in-in path integral formalism. This approach shows clearly that 
the tunneling process becomes strongly coupled whenever the sound speed of the tunneling 
trajectory departs sufficiently from unity. We argue that these two very different perspec- 
tives are nevertheless consistent, and may imply the existence of a simple resummation of 
the strongly coupled interactions of the field. 
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1. Introduction 

Tunneling is the best understood non-perturbative phenomenon in the theory of gravity 
coupled to scalar fields. It is central to the description of vacuum phase transitions, dy- 
namics on the landscape, and the initial conditions for inflation. There are two types of 
tunneling in nature: quantum tunneling and thermal tunneling. Quantum tunneling is a 
classically forbidden process, while thermal tunneling appears for any open system whose 
environment induces a finite temperature and thus can be purely classical. In any finite- 
temperature tunneling system, these processes coexist. The precise nature of the potential 
barrier through or over which a system is tunneling and the temperature prevalent in the 
system determine which kind of tunneling is dominant. Thermal tunneling dominates at 
high temperatures and for very broad potential energy barriers, whereas quantum tunnel- 
ing dominates for narrow barriers. When considering tunneling in cosmological theories 
with a positive potential, the effect of gravity on tunneling is equivalent to adding tem- 
perature. This is because of the finite temperature seen by an observer in de Sitter space. 
The presence of gravity thus implies the coexistence of quantum and thermal tunneling in 
cosmology. 

Many of the novel inflationary models proposed in the past few years consist of a 
scalar field with a non-minimal kinetic term - at least at the level of their 4-D effective 
field theory. The example with the best physical motivation is Dirac-Born-Infeld (DBI) 
inflation [p. In this model, the inflaton is the scalar field describing the location of a 
D-brane in a warped throat, and the action is of a square root form characteristic of 
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Figure 1: A typical tunneling potential. The minimum labeled (/)taisc denotes the false vacuum, 
while (/)truc is the true vacuum, separated by a barrier with maximum at ^barrier- 

relativistic particle motion. More generally, A;-inflation and /c-essence are examples 
of field theories of this type. In such theories, the Lagrangian is an arbitrary function 
of the scalar field and the usual kinetic term £ = yJ—gp{X, (j)) with X = —^(dcf))'^. In 
fact, all possible predictions of single field inflationary models coming from local theories 
can be described within a unified framework using this kind of approach [Q]. In any case, 
non-minimal kinetic term models represent explicit, computable alternatives to standard 
slow roll inflation. 

In the case of DBI inflation, quantum tunneling - i.e., tunneling within the Coleman- 
de Luccia framework - has been considered in Ref. In this work we shall focus our 
attention exclusively on thermal tunneling. As we discuss below, there are two ways to 
describe thermal tunneling: a local description and a global description. Which one we 
choose depends on whether we wish to study the physics seen by a single observer, or want 
to ask about the global geometry, including superhorizon modes. The local description 
is straightforward to generalize to a variety of models, but generalizing the global picture 
proves to be a more challenging undertaking. There is an immediate, technical challenge: 
the slow roll approximation cannot be used. This limitation requires a generalization of 
the stochastic description of inflationary fields; this we carried out in a previous paper 
Q and was addressed in earlier work 0, |8|. There is, though, a more fundamental hur- 
dle: the presence of intrinsic non-Gaussianity in these theories plays a crucial role in the 
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tunneling process. This may prevent any perturbative calculation from capturing the full 
phenomenon, forcing us to develop new techniques. 

1.1 Global versus local 

There are two ways to interpret thermal, also known as stochastic, tunneling. To un- 
derstand these two approaches most clearly, it is helpful to work in the decoupling limit 
(Mpi -^> oo) with fixed Hubble parameter, H. In this limit, the scalar fields are taken 
to exist on a fixed de Sitter background geometry; that is, gravity does not respond to 
the changes in the scalar fields. We begin by considering the local approach. The local 
approach describes the tunneling process from the perspective of an observer within a par- 
ticular causal patch. This patch is covered by what are known as the static coordinates in 
de Sitter: 

ds"^ = -(1 - H^r^)dt^ + (1 - H^r^y^dr^ + r^d^n. (1.1) 
In this local, or subhorizon, approach, the natural coarse graining of the closed system is: 

• Open System: causal patch, sub-horizon modes, defined in position space as fields 
with support in r < H^^ 

• Environment: super-horizon modes, defined in position space as fields with support 
in r > H^^ 

• Tunneling Rate: Determines rate that the observer's Hubble-sized region "jumps" 
over the barrier 

The only effect of the environment is to induce a finite temperature T = H/2tt for the 
subhorizon system Q. This idea is known as de Sitter complementarity |jl^. Consequently, 
in this picture tunneling can be described by a finite temperature field theory living in the 
causal patch. 

In the global approach, we consider physics from the perspective of a 'super observer' 
who can see physics on a global slicing. For simplicity, we can take this to be the flat slicing 
in de Sitter, 

In the global, or superhorizon, approach, the natural coarse graining of the closed system 
is: 

• Open System: Long wavelength modes, defined in momentum space as fields with 
support on length scales A > e^^aH~^ ^ 

• Environment: Subhorizon modes, defined in momentum space as fields with sup- 
port from length scales A < e~^aH~^ 

• Tunneling Rate: Determines rate that the field, coarse grained over several Hubble 
volumes in flat slicing, "jumps" over barrier 

^Here e is a small parameter introduced so that the the coarse graining scale is a few Hubble radii where 
the growing mode has already come to dominate over the decaying mode. 
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The effect of the environment here is to induce a stochastic noise for the long wavelength 
system. Adding gravitational dynamics, whilst complicating the description of the system, 
does not significantly change these considerations. 

We stress that these two methods are not calculating precisely the same tunneling 
rate. Nevertheless, they are sufficiently similar that we expect a comparable result. In 
brief, they both answer the question: 'What is the rate at which we should expect to find a 
Hubble volume exiting a metastable vacuum and reaching a stable - or at least more stable 
- vacuum through a stochastic/thermal fluctuation?' 

2. The local picture: Euclidean Quantum Gravity 

In the local picture, we consider the causal patch to be at finite temperature T = H/2'k. 
The system is thus described by a partition function 

Z = Trs e"^^='='«^ (2.1) 

with (3 = 2itH~^. The justification for this is as follows. If the global vacuum is chosen 
to be the Bunch-Davies vacuum |0), then the state that an observer sees is obtained by 
taking the trace of the density matrix over those modes with support outside the causal 
patch. This is is found to be equivalent to a thermal state [^: 

PS = T¥l (|0)(0|) = z-^e-^-^^'-'-. (2.2) 

Here -ffstatic('?i"; (p) is the static patch Hamiltonian (written as a function of the field and its 
spatial derivatives as well as the conjugate momentum) which is the generator associated 
with the timelike Killing vector in de Sitter. Standard techniques tell us that evaluat- 
ing this partition function is equivalent to calculating the Euclidean path integral over 
configurations which are periodic in Euclidean time te with period /?. That is, we may 
write 

Z= D^Dtt e-^i^'^'^'''^ (2.3) 

where ^ 

5^(0, vr) = y dTE y drj r^d^n (^vr^ - Wstatic(vr, 0)) . (2.4) 

Note that this is a canonical representation of the path integral, which is necessary to get 
the correct measure in the case of a non-minimally coupled field. For a minimally coupled 
scalar field, the explicit form of the static Hamiltonian is 

^static = y dr j r^d^n ^static (2.5) 

= Jj^drj d'n Q(l - H\') + r'^l^ + i(Vn<A)' + rV(<^) 
For a configuration in which (p is fixed at a minimum or a maximum, the energy is 

i^static = ^^Vi(l)). (2.6) 
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When the field has reached a true equiUbrium, the above noted partition function is reah 
tunneling from the left of a barrier to the right of a barrier is exactly balanced by the 
reverse process. Nevertheless, we can infer the rate for tunneling from the false to the 
true vacuum by restricting the range of the phase space path integral so that it only 
includes configurations with right moving flux and in which the initial state is prepared 
as being thermal to the left of the barrier and vacuum to the right. This path integral 
can be evaluated in the saddle-point approximation. Tunneling instantons are those whose 
fluctuations exhibit an odd number of negative modes. To deal with these, we must perform 
a Wick rotation of the integration contours. This Wick rotation picks up factors of i, and 
the total tunneling rate is determined by the imaginary part of the partition function 

rocJmlnZ. (2.7) 

A class of instantons that always exist are defined by the static configurations 

(^i^static(0,7r = 0) 

TT = 0, — = 0, (2.8) 

where the field is fixed at a maximum or minimum of its potential. Unlike in Minkowski 
spacetime, these are guaranteed to have finite action in de Sitter spacetime because the 
static patch has finite volume. When the instanton is at a maximum, there is at least 
one negative mode, hence there is a nonzero contribution to the tunneling rate. In the 
decoupling limit, the tunneling rate is given by the Arrhenius law 

r OC exp {-13 (-frstatic(0barrier, 0) - i?static('/'false, 0))) (2.9) 

where (pMse is the false vacuum and (/^barrier is the maximum of the potential intermediate 
between the false and true vacuua. This is precisely the decoupling limit of the Hawking- 
Moss result 

r CX exp (S'cnt(0barrier) " Sent (^false) ) , (2-10) 

where Sent is the entropy of the de Sitter geometry associated with the extrema: 

S,M) = ^ = (2.11) 

This follows since 



A5ent ^7^^^ ~ -27TH-' I-——AV ] = -/JAi^static (2.12) 

In a similar way the Coleman-de Luccia instanton can also be interpreted as finite temper- 



ature quantum tunneling |11] in the decoupling limit. 



2.1 Extending to non-minimal models 

How does the local picture change in a theory with a non-canonical relationship between 
momentum and kinetic energy in its action? Let us consider theories of the general form, 

S^= I d^x^gp{X,<t>). X = -\{d4>f (2.13) 
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The best known example of such a model is DBI inflation. There, the pressure p takes the 
explicit form 

cP) = -T(0) - 2X/T{cP) + T{4>) - (2.14) 

where T[(j)) is the tension of the D-brane moving in the warped throat. The chief novel 
physics encoded in this action are that perturbations travel with a sound speed less than 
the speed of light. This sound speed is given by 

= 1^ = (2.15) 

dp P,x +2Xp,xx 

(note that here and throughout the paper we will use commas to indicate partial differen- 
tiation; that is, P,x = dP/dX). This means that when we perturb around a cosmological 
solution, the perturbations travel at a speed typically lower than unity, and there is a sound 
horizon which is smaller than the causal horizon. In the local picture we are discussing, 
this reduced sound speed implies that Lorentz invariance is broken for the scalar field vacua 
described by such theories all the way up to the cut-off energy scale when this effective 
Lagrangian is expected to breakdown. To see this, let us consider the short distance singu- 
larity structure of the Feynman propagator for fluctuations 5(j) around some cosmological 
solution, taken to be in their adiabatic vacuum. Since this will be dominated by the con- 
tribution from subhorizon modes, it is sufficient to treat Cg as a constant, since it varies 
'slowly'. Nevertheless, even high momenta modes travel at the sound speed Cs, and so the 
Feynman propagator will have the Lorentz-violating form 

hm (r50(x)<50(x')) = L 1 7 =2 ..^ ^,.o . ^ (2-16) 



Svr^CsP,^ J {t - t'Y - Cs"^ ex.p{2HT) {x - x'Y - ie 
as if the scalar field perturbations were living on the emergent geometry 

ds^ = {csP,x) {-dT^ + cj^e^^^dx^) . (2.17) 

In the case where Cg and p^x are exactly constant, the emergent geometry is also de Sitter, 
with Hubble constant H/ ^CsP,x and associated emergent background temperature equal 
to Tg = H/ {27r^CsP^x)- In the case of a DBI-like kinetic term p^x = and so there is 
no change in the apparent temperature within the emergent geometry. The horizon of the 
emergent geometry is the sound horizon which the perturbations see and is in general no 
longer the same as the true event horizon. 

Nevertheless, the physical temperature of the background is still the same. Following 
the argument of Ref. j^, if the Feynman propagator is defined to be analytic in the emergent 
geometry for spacelike separated points then it will be invariant under 

T^r + t-, (2.18) 

which is the same as requiring analyticity of the propagator in the real geometry. This 
is the key step in justifying the notion that subhorizon physics in de Sitter is equivalent 
to a finite temperature system of temperature T = H/2'k. This is true when Cg and p^x 
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are constant because the emergent metric depends on time only through e , which is 
invariant under this shift. 

What this example illustrates is that even if we choose the global state to be the adia- 
batic vacuum state for fluctuations around nontrivial cosmological solution, with nontrivial 
speed of sound, i.e. even if we choose the global vacuum not to be Bunch-Davies, it still 
makes sense to think of the sub horizon system as being at a finite temperature T = H/2'k 
set by the event horizon and not the sound horizon. 

In short, causality demands that the part of the global vacuum describing superhorizon 
physics is never relevant to a local observer. Let us denote the Bunch-Davies vacuum by 
|0). At a given time, the apparent vacuum la) for perturbations around some cosmological 
solution formally corresponds to a coherent state^ 

|a) = U\{)) = e-l"l' exp[aa + a^a^] |0), (2.19) 

where the a and are defined for the particular state, and aa is shorthand for a sum 
over a complete set of modes ^^ajOj. We can always divide the creation operators into 
superhorizon (L) and subhorizon (S) modes, 

a = aL + as, (2.20) 

such that we guarantee [01,03] = [ol, a|] =0. As a density matrix this state is 

\a){a\ = C/|0)(0|C/-i = ULUs\Q)ms^Ul\ (2.21) 

where Us = e"!"^!^ exp[a5a5 + 0:505] and similarly for Ul- ^ local observer can only see 
the part of the state obtained by tracing out the long wavelength modes, namely 

PS = TVl|«)(«| = Tii^iULUsmms^Ul^) (2-22) 
= Us{'Tii^UL\Q)m^^)Us^ = C/5(TrL|0)(0|)t/5i (2.23) 

= UsZ-'^e'l^^^^^^^'^Ug'' = Z-^e"^'^«-'^«*^«'=^s\ (2.24) 

In other words, if the global vacuum is the adiabatic vacuum defined around some cos- 
mological solution, the state seen by a local observer is a coherent excitation of a thermal 
state with the same temperature T = H/2tt. 

Nothing about these changes implies that the usual thermal tunneling instanton should 
no longer exist or be substantially changed. By construction, the thermal instanton in this 
picture occurs via a static field configuration. Since the novel effects from the nonminimal 
kinetic term only arise when velocity is finite (for instance, in DEI inflation, p{X, (j)) ~ 
X — 2X^/T{(f)) — V{(f)) for (p <C VT), there is little influence of the kinetic term on the 
static configuration. Thus, the instanton itself and the fluctuations around it - which are 

^Technically speaking, the coherent state may be in a unitarily inequivalent representation since the 
sum |Qip may not converge in the UV. In the context of field theory with a cutofT, this distinction is 
irrelevant, and is really a reflection of the fact that we are not allowed to excite modes lying beyond the 
cutoff. We take the point of view that these manipulations are meaningful if, having defined them for fixed 
cutoff, they make sense in the limit that the cutoff tends to infinity, after appropriate renormalizations. 
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necessary to determine the existence of a negative mode and hence the existence of the 
instanton - are not affected by the nonminimal kinetic terms. 

Any difference between the nonminimal and minimal cases will come from higher-order 
corrections, specifically two-loop corrections in the Euclidean path integral. The instanton 
calculation is a saddle-point calculation, usually considered to one-loop order, that is only 
valid as long as the two-loop and higher corrections are negligible. Fluctuations around 
the instanton generate terms in the Lagrangian schematically of the form 

5C = p,x6X + ^p,xxSX^ + ^p,<f,^<p5<P^ + p,x^5X6(t> = ... (2.25) 



where 5X = -\{d5(t)f . A necessary, but not sufficient, condition for the one-loop calcula- 

5X < 



2 

tion to be valid is 



P,xx 



(2.26) 

barrier !^ 



since at the point of equality the fluctuations in the higher order contributions to the 
kinetic term become as large as the quadratic part. For a minimally coupled scalar field 
p^xx = and so this condition is always satisfied. In the language of instantons, when 
this condition is not satisfied, the fluctuation determinant can no longer be determined by 
a gaussian integral about the instanton solution. 

The typical magnitude of the fluctuations is set by the Hubble scale 5(j) ^ H , dcp ^ H^. 
In the case of a DBI model, this condition for the breakdown of the one-loop approximation 
is, then, 

~ — ~ ~ 1. (2.27) 

-^('/'barrier) ^(^barrier) ^(^barrier) 

This result is intuitive. We expect loop corrections to the effective action to be of order 
at least H^. On the other hand T{(l)) sets the scale of the kinetic term for DBI, so this 
is just the condition for the loop corrections to the effective action to become comparable 
to the classical kinetic term. This is, of course, not the only condition for the validity 
of the one-loop approximation. However, the additional conditions have to do with the 
(j) dependence p{X, (j)) and so are typically easier to satisfy if the potential is sufficiently 
broad. 

This is similar but not the same as the perturbative bounds that have been placed 
on eternal inflation, when quantum fluctuations dominate over classical evolution in the 
field evolution in Refs. |]l^, The key difference is that in our case we are expanding 
around a static solution (the instanton), whereas in these other works the authors have 
found nontrivial bounds by expanding around a time dependent solution with Cg ^ 1. The 
two calculations are qualitatively similar in that they reflect the break down of pertur- 
bation theory for the fluctuations - in one case around the instanton, in the other the 
cosmologically rolling solution - but should not be expected to match in detail. 

We have thus established a subset of cases where the thermal tunneling instanton is 
not valid. It will be a bad approximation when the top of the potential barrier is at such 
a high energy that all perturbations are dominated by relativistic effects. In contrast, the 
Hawking- Moss instanton for a minimally coupled scalar field is valid |15] even when the top 
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of the potential barrier reaches the energies associated with eternal inflation. This is closely 
connected with the observation that a minimally coupled field reaches eternal inflation 
within the regime of perturbative control ||l^, but a non- minimal theory such as DBI 
inflation is not under perturbative control in the regime of eternal inflation, at least when 
the sound speed departs sufficiently from unity. The difference is that the loop expansion 
for a minimally coupled field is determined by the potential in the decoupling limit, rather 
than by its kinetic term. If the potential is sufficiently flat, then these corrections are 
always negligible. This is always true when slow roll is a good approximation, and can be 
true even at the top of a broad potential barrier in a tunneling set up. 

We now turn to the complementary, global picture of de Sitter space and inflation, 
where the techniques of stochastic inflation can help us to understand better what exactly 
is going on in thermal tunneling. 

3. The global picture: stochastic inflation 

In any accelerating geometry, perturbations of a quantum field expand faster than the 
Hubble radius. As these perturbative modes cross the Hubble horizon, their two quantum 
oscillating components split into a growing and a decaying classical mode. The long wave- 
length state of the quantum field retains any non-linearities it possesses, but is squeezed 
by the expansion and becomes progressively more classical as its growing mode dominates 
over its decaying mode. The long wavelength or background field evolution thus has two 
contributions: part from its overall classical evolution, part from the stochastic kicks it 
receives as each subhorizon mode exits the horizon and adds itself to the collection of 
superhorizon modes constituting the background. 

This process is always operative for infiating geometries, but the stochastic influence of 
the subhorizon noise is generally small enough to ignore. However, the process of thermal 
tunneling is precisely the exception that proves this rule: a stochastic jump of rare vigor 
comes along that is of sufficient size to move the average field at its false vacuum value in 
one Hubble patch over the potential barrier that separates it from the true vacuum. The 
other time when large perturbations are important is, of course, during eternal infiation. 
Both of these phenomena can be described within the framework of stochastic inflation. 

3.1 Deriving the stochastic formalism 

A rigorous derivation of the stochastic approach is quite involved, but we can sketch the 
details as follows: Consider the evolution of the system in the flat slicing with canonical 
action 

5[P, </.] = y drd'x e^^^ (^P^ - Wflat(P, <A)) , (3.1) 

The momentum conjugate P defined here should not be confused with vr which was the 
momentum in the static slicing. For instance, for the case of a massive minimally coupled 
scalar field on de Sitter we have 

Wflat(P, <P) = \P^ + ^(V<A)' + \m^\ (3.2) 
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with a(r) = e^"^. Expectation values of fields can be derived from considering the in-in or 
closed-time-path integral, which contains two branches + and — . For instance we have 

(0|Te*^-^'^|0) = j D(i)+DP+D(i)^DP^e'^^P+^^+^-'^^^'''i'~^+'^'-^^+, (3.3) 

and the path integral is performed with final time boundary conditions 0+(rj) = 0_(rj), 

P+{Tf)=P.{Tf). 

To connect with the usual stochastic formulation, we want to construct a probability 
distribution on phase space. In an appropriate limit, this approach will give rise to a 
Fokker-Planck equation describing the stochastic evolution. In quantum theory, for a 
density matrix p we can define a quasi-probability distribution - the Wigner function - as 

follows /• , , f f 

W\4>,P;t\ = J D£e'f'''--'"^P^{cP--\p{TM + -). (3.4) 

As is well known, this function can be used to calculate expectation values in the normal way 
as if it were a probability distribution defined on phase space. The quantum information 
is contained in the fact that W is not positive definite. 

The Wigner function can also be thought of as the 'wavefunction' associated with the 
functional Fourier transform of the path integral ( |3.3| ), given the replacements (p± = (p^£/2 
and P± = P± x/2 

W[(t>,P] = j"^ Vcj) j VP j V£ j Vxexp(^iS^ , (3.5) 

with 

S = f drSxe^^^ 



Xdr(t> + £ {drP + SHP) - Wflat[(/' " |' ^ + f ] 
+Wflat[0+|,P-|] 



(3.6) 



Just as in the local approach, we traced out the superhorizon degrees of freedom to deter- 
mine the subhorizon density matrix, here we define the long wavelength density matrix 
by separating each field A = (</), P, x) into long and short wavelength modes 

Aix)=ALix)+As{x), (3.7) 
-'^■-<P{k,T)W{k-T), As= -—:.e''-^cl>{k,T){l-W{k;T)). (3.8) 



Here W{k;T) is the Window function that distinguishes super and sub-Hubble modes, and 
is most commonly taken to be of the form 

W{k;T) = e{eaH -\k\), (3.9) 

with e a constant whose value is taken to be much less than unity. This is in fact a 
crucial point: the coarse graining is not performed at the Hubble horizon, because the 
state still appears quantum there. The state only appears classical once the squeezing 
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mechanism has taken effect, when the modes are a few e- folds outside the Hubble horizon. 
This split is defined so that equal time bilinears of the form (for some functions A and B) 
J d^xAi^(x, t)Bs{x, t) = 0. However since the coarse graining is explicitly time dependent, 
we have 

J d^xALdrBs = - j d^x{drAL)Bs = - J -^A{-k)B{k)W{k;T)drW{k;T) 
ff] H"^ f d^k 

= 2~ J J^Ai-k)B{k)6i\k\ - aH). (3.10) 

This means that independent of the form of the Hamiltonian, there is always a direct 
coupling between the momentum and field of the short and long wavelength modes through 
the Pdr(f> term in the Lagrangian. For a free field theory, this is the only coupling term. 

We construct the reduced Wigncr function describing the effective long wavelength 
dynamics by tracing out the short wavelength modes 

Wr [4>L, Pl] = J V(t>s VPs W [0i + 4>s. Pl + Ps] ■ (3.11) 

The reduced Wigner function is the 'wavefunction' associated with the path integral, 

Wr [ct>L, Pl] = j Vct^L VPl V8l Vxl exp (iS [cf)L, Pl, £l, Xl]) F [c^l, Pl, £l,Xl] (3.12) 

where F is the Feynman- Vernon infuence functional 

F [(f>L, Pl,£l, xl] = J 'D(f>s VPs V£s Vxs exp (iS [<f>, P, £, x] - iS [(f>L, Pl,£l, Xl]) ■ 

(3.13) 

The full expression for F is obviously difficult to determine, however in most cases of 
physical interest, logF has a controlled perturbative expansion in powers of £l and XL- 
We will not require this sort of expansion for our calculations. However, let us write down 
concretely what form such an expansion would take. Since we have F* [(f)L, Pl,£l,Xl] = 
F [(pL, Pl,£l, — Xl] , we can simply expand log F in the following way: 

log F = C [<i>L, Pl] + j {iA^ [H, Pl] £l + iAp Pl] Xs) e^^'^d^xdr (3.14) 

-JJ a!Wye3^(^-+^^) Q^l(x)D<^^(x, y)^L(2/) + ^^l(x)D<^p(x, y)xL(2/)+ 

Ixl{x)Dpp{x, y)xL{y)^ + 0{£^ x\ £\, £x^) , 

For some functions C, Ax, and Dxx that will encode the particularities of the system 
one is studying (where X stands for the field or momentum, as above) . The influence 
functional can also be represented by its functional Fourier transform 

F [H, Pl, £l, xl] = j Vn^VripV [r?^, vp, H, Pl] e-^/K^^+^^^^]-'^^<^'-'^- , (3.15) 
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with and rjp being the conjugate variables to £l and xl- We will see soon that these 
r]x will function as noise forces in the equations of motion for our system. The advantage 
of this approach can now be seen. In the semi-classical limit we can approximate 



S[(^l,Pl,£l,Xl] ~ J d^xdre^''^ 



XLdrC^L + £l [dr + 3H] Pl + " 



where Hl = f d^x7iaa.t{4'L, Pl)- We can then easily perform the path integral (Eqn. |3.12| ) 
over £l and xl with the result 



r [(PL 



,Pl] 



'>L fPL 

V(j)L / -DPl I Vj]^ I Vr]pV[n<t>,VP,<pL,PL] 



X 6 



drC 



6Hl 
6P 



drPL + 3HPl + ^ - r?p 



(3.16) 



Thus we see that the reduced Wigner function is the classical probability distribution 
associated with the two-noise Langevin equations, as considered in reference @ 



dT4>L 

drPl 



6P 



+ V<t> 



-3HPl J— + VP 



(3.17) 
(3.18) 



where V in Eqn. 3.16 is a quasi-probability distribution that encodes the "coloring" - the 
statistical behavior - of the noise terms. 



3.2 Irreversibility and the fluctuation theorem 

Having connected with the approach of we can now address the question of how to 
calculate the tunnelling rate. The key observation is that the dominant contribution to the 
tunneling rate is determined by the irreversible nature of the coarse-grained dynamics. A 
crucial difference between the stochastic approach and the local approach of the previous 
section is that in the stochastic case the effective dynamics of the long wavelength system 
is irreversible due to the presence of dissipation (Hubble damping) and diffusion (noise 
from subhorizon quantum fluctuations). In the local case causality guarantees that the 
superhorizon and subhorizon modes decouple, and hence there is no diffusion, and since 
the static time coordinate is a Killing vector, there is also no damping seen by a static 
observer. 

To be more precise about what we mean by irreversibility, consider the propagator /C 
that evolves the reduced Wigner function between two times Tj and tj 



Wr 



Vcf>lvPlic cf>{,p[,Tf,<pi,Pln 



.Pin]. 



The propagator /C has the path integral representation 

lc[<p{,P[,Tf;(^i,Pi,Ti = 
rPL{Tf)=P[ 



(3.19) 



(3.20) 



r<t>L{rf)=cl>,^ 

I Vc^ j VP I V£ I VxL e^p(iS[^L,PL,£L,XL]) F [4>l,Pl,£l,Xl] 

4>L{Ti)=4>l JPL{n)=Pl 
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The fundamental statement of the irreversibihty of the system can be encoded in the 
following ratio 



7^ 



P[,Tf,<pi,Pln : (3.21) 



JC 



In other words, for a given transition between two states vl — > i?, we can compare the 
probability for this transition, with the probability for the time reversed transition B ^ A 
to occur in the same time. For a reversible system, we have 7^ = 1, and so TZ encodes any 
departure from reversibility, which in practice is a measure of the amount of dissipation 
and diffusion. 

A priori, this ratio can be a complicated function of the initial and final phase space 
values and times. However, in many well understood systems, the ratio takes the simple 
detailed balance form 

n = exp (m[c/>{, p/] - MWl, Pi]) , (3.22) 

for some function, M. This result is closely connected to a result known in the condensed 
matter literature as the fluctuation theorem [16|. It states that the ratio of the probabilities 



of making a transition from A ^ B to the ratio for the time reversed process — > j4, in 
some time r/ — Tj is the exponential of the amount of entropy produced in the transition: 

where Sent is the entropy associated with the system. Thus it is natural to identify 
M[(/>^,P2] ^ the entropy associated with a given phase space configuration. 

The fluctuation theorem is the modern statement of the Second Law of Thermodynam- 
ics, encoding the fact that the probability for entropy to increase is greater than that for it to 
decrease, recognizing that whilst spontaneous entropy decrease is allowed, its probability is 
highly suppressed. This exponential suppression for the reverse process, directly determines 
the exponential suppression of the stochastic tunneling rate. The rate to go from the false 
vacuum to the top of the barrier is given by the exponential of the de Sitter entropy differ- 
ence between the two configurations multiplied by the rate to go from the top of the barrier 
to the false vacuum: P(false — > barrier) = exp(S'ent(<^top) — 5'ent(0faise))-P(barrier — > false). 
The latter process P(barrier false) though, is a classically allowed process; hence its 
transition probability will contain no additional exponential suppression. Similarly, once 
the field has reached the top of the barrier it can roll classically down towards the true 
vacuum. Thus we infer that the tunneling rate is given by 

r(false true) OC exp (S'ent ('/'barrier) " Sent{<P{a.lse)) , (3.24) 

up to factors determined by the classically allowed processes. 

So, how does this connect with cosmology? We have performed the above analysis for 
inflation with a minimally coupled scalar field. We have expanded the influence functional 



only to quadratic order as in equation |3.14|. In this case, one can explicitly show that in 
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the ultralocal limit, the detailed balance condition [ 3.22 | holds, and furthermore that it is 
consistent with the fluctuation theorem: the change in M is precisely equal to the decou- 
pling limit of the change in the entropy (eq. p.l2|| ) of a scalar field in a de Sitter spacetime 
determined by the potential V{(j)). Even when gravitational backrection is included, this re- 
sult holds up |6|, |l8|. In other words, the dominant contribution to the stochastic tunneling 
rate is identical to the rate determined by Hawking and Moss in the local case via the in- 
stanton method. This strongly supports the interpretation of the Hawking-Moss transition 
as being equivalent to the stochastic tunneling process, as is usually assumed [p^] . 

Our interpretation of this result is that the two results agree in the usual case (Hawking- 
Moss and stochastic) precisely because of thermodynamic reasons. In the Hawking-Moss 
case this is explicit, because it is a thermal calculation of the imaginary part of the free 
energy. In the stochastic case, it is explicit because the evolution of the Wigner function 
is consistent with detailed balance and the fluctuation theorem. 

3.3 Irreversibility for DBI inflation, non-gaussian corrections to tunneling rate 

In a previous paper |^] we extended the stochastic formalism to DBI inflation, and in 
the language of this section, derived the Langevin/Fokker-Planck equations that follow 
expanding the log of the influence functional to quadratic order as described in equation 
| 3.1^ , and then taking the ultralocal long wavelength limit. We found that to this order, 



the evolution of the (Wigner) probability distribution on phase space, indeed satisfied 
detailed balance, in the form 

n = exp (m[</){, p/] - M[0i, Pi]) , (3.25) 

but that the full expression for M was 

M[0,P] = M[<A, J] = -IGvr^l d<kMl, (^^^^) (3-26) 

where J is related to and P via the canonical transformation used to derive the Hamilton- 
Jacobi formalism Q. As pointed out in Refs. |6|, |l^ if this result is taken at face value, 
the presence of the additional factor of Cg inside the integrand means that the associated 
stochastic tunneling rate, whose suppression is by the above arguments determined by IZ^ 
is no longer given by the same result as the local thermal tunneling rate, since M is no 
longer identical to the expression for the entropy of de Sitter spacetime, with the associated 
Hubble constant H{(j),J). It is only in the case Cs = 1, in other words in the minimally 
coupled scalar case do we recover the usual answer. 

The stochastic tunneling rate and the Hawking-Moss tunneling rate for models in which 
Cg may depart from unity along the tunneling path thus appear to disagree. One might 
think that this apparent discrepancy should be taken as evidence that the two rates do 
not describe quite the same thing. We find that result unlikely since, as we explained at 
length in the introduction, they both effectively describe the rate for a Hubble sized region 
to traverse the barrier between the false and true vacuum. If this is not the case, then 
it must be that the approximations used in deriving one of the two results breaks down. 



- 14 - 



In the local instanton method, we have already argued that the calculation is usually well 
under control from higher loop corrections to the tunneling rate, as long as the condition 
^ r((/)barrier) is Satisfied. This suggests that it is the stochastic calculation that is 
failing. 

Let is consider the approximations used in deriving the stochastic formalism, sketched 
out above. The key approximations are: 

• Long wavelength perturbations are treated as classical: This is encoded in equation 



[3.16| where S is expanded to linear order in £l and xl- 



Long wavelength modes decouple, i.e. the ultralocal limit: In this approximation we 
can neglect spatial gradients for the long wavelength fluctuations, and the noise 
correlators are assumed to be local in space (r/(x)r/(x')) oc 5^{x—x'). This is explained 
in more detail in Ref. [01. 



• Short wavelength noise can he treated as gaussian: This is encoded in equation 3.14 | 
where the influence functional is taken to be gaussian in £l and xl-) 

The first and second properties follow from the squeezing of the long wavelength quan- 
tum state and causality respectively. Both of these properties are equally well satisfied in 
DBI inflation as they are in usual slow roll inflation. However, the third approximation is 
not, precisely because in DBI inflation, the quantum fluctuations in the inflaton are more 
non-gaussian, then in turn the stochastic noise which arises from the short wavelength 
fluctuations is more non-gaussian. 

Our central claim is that it is the neglect of the non-gaussian noise fluctuations that 
is responsible for the apparent discrepancy between the stochastic tunneling rate and the 
Hawking- Moss result. As soon as Cg departs noticeably from unity along the tunneling 
trajectory, the non-gaussian contributions to the irreversibility ratio IZ become important. 
Indeed it seems plausible that the physical consistency of these two results requires that 
the ratio IZ is always consistent with the fluctuation theorem, and so when the corrections 
from the non-gaussian noise fluctuations are resummed, we will recover the usual result. 

To restate this problem for DBI inflation, whenever Cs departs from unity along the 
tunneling trajectory, perturbation theory in the variables £l and xl which are conjugate to 
the noise in the in-in path integral breaks down. The path integration over these variables 
must be performed non-perturbatively. Unfortunately it is not clear how to achieve this as 
yet. 

To see that the noise fluctuations are becoming non-perturbative during tunneling, let 
us compute the three-point function for the noise ry^. Following the coarse-graining recipe 
described above, and in more detail in Ref. P, the noise ry^ is related to the subhorizon 
fluctuations as follows (note: in the ultralocal limit we can drop the x dependence) 

ri^ir) = - j -^,a{T)HH{\k\ - eaH)Mk). (3.27) 

where (j)s are the short wavelength fluctuations about the long wavelength background (pi- 
The key point is that the noise has a non-zero three-point correlation function, precisely 
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because the subhorizon fluctuations have a non-zero three-point function. As we are only 
interested in the order of magnitude, we can estimate the three-point function as follows: in 
the framework of cosmological perturbation theory we would define the comoving curvature 
perturbation ( ~ H(f)s/4>L- It is usual to encode the non-gaussianity of the the fluctuations 
in terms of the parameter /jvl which is in general scale and shape dependent, 

{C')--lfNL{{Of. (3.28) 
In terms of the fluctuations <ps this is 

(4)~-^/iVL((05))'. (3.29) 

(PL 

Given the expression for the noise in terms of the fluctuations (ps, we find the three-point 
function for the noise to be of the form 

{V4>{n)V<t>{T2)V4>{'r3)) ~ fNL^S{Ti - T2)S{ti - T3), (3.30) 

<PL 

where we have made the simplifying assumption that the mixing of modes that exit the 
horizon at different times is negligible, which provides the two delta functions. This should 
be compared with the form for the two-point function for the noise 



(??0(n)r?^(T2)) - i^^5(Ti - T2). (3.31) 



Consider the dimensionless noise variable o(r) = ^ which encodes the ratio of the stochas- 
tic fluctuations to the classical evolution. We can average this over the tunneling trajectory 
that takes time At 



At-' 

lo 



/ dTq{T) (3.32) 
Jo 

then a simple order of magnitude estimate gives 



~ S^At (3.33) 

and 

{q') - {fNL)S^Ar' (3.34) 

Here we understand {/nl) to be the average of /nl along the tunneling trajectory. In other 
words just as {/nl) naturally encodes the nonlinear ities in it also naturally encodes the 
non-linearities in q 

{f)-{fNL){{f))'. (3.35) 

Now the key point is that since q measures the ratio of the stochastic fluctuations to the 
classical evolution, during tunneling, which is a rare quantum event where the quantum 
fluctuations are sufficient to dominate over the classical evolution, wc inevitably have g ~ 1. 
But unless {/nl) ^ 1, for any such large fiuctuation all the higher order correlation 
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functions of q become of order unity, and hence perturbation theory breaks down in the 
noise variable. 

We have expUcitly checked by a more detailed calculation that the above argument 
goes through when the corrections to the irreversibility ratio IZ are computed from the non- 
gaussian noise fluctuations. The bound on the perturbative region we find for tunneling, 
namely {/nl) Ij is much more restrictive than that discussed in Refs. 0] for 
vacuum fluctuations about a cosmological solution. This is not surprising since we are 
looking at different effects. Tunneling is itself a non-perturbative phenomena, and so 
the perturbative bounds derived there are not appropriate. Tunneling is described by 
rare stochastic fluctuations, and the conditions for the validity of the stochastic tunneling 
calculation are necessarily more restrictive for these rare fluctuations, in comparison to the 



typical fluctuations about cosmological solutions considered in Refs. [13, 14]. 

Returning to the tunneling rate, which as we have argued is principally determined by 
the irreversibility ratio TZ we find that the correction to InF from the non-Gaussian term 
takes the form 

51nr~ (/jvL)lnr, (3.36) 

with similar contributions from higher A^-point functions. To get a reliable result, we 
would have to sum this full set of contributions. For minimally coupled scalars we always 
have {/nl) ~ (fisiowrou) 1 and so these higher order corrections are negligible. This 
explains why in the standard case there is perfect agreement between the local and global 
descriptions of the tunneling process ||l8[| . 

3.4 A contradiction? 

The fact that the effect of the non-gaussian noise correlation functions become significant 
during tunneling is evidence that the dynamics of tunneling in DBI inflation is strongly 
coupled, i.e. strongly sensitive to the higher order kinetic terms. Since the strongly coupled 
regime is difficult or impossible to treat perturbatively, we should not be surprised that a 
perturbative expansion of the logarithm of the influence functional is insufficient to capture 
these dynamics fully. In practice this means that the familiar picture of stochastic inflation 
with gaussian noise is no longer appropriate for DBI inflation, at least in the description 
of tunneling. 

However, the Euclidean calculation gives the tunneling rate independent of the tun- 
neling trajectory's dynamics to be 

r ~ exp (S'cnt(0barrier) — Sent ('/'false)) , ^ent = ~Jp(^ ' (3.37) 

This result is identical to that for inflation driven by a minimally coupled scalar field and 
seems to be insensitive to the nontrivial kinetic terms in the DBI action, at least as long 
as we satisfy H'^ <C T(</)barrier)- If these two calculations are really describing the same 
physics, it is difficult to see how it could be the case that one rate is essentially insensitive 
to the non-minimal kinetic terms, whilst the other calculation is strongly sensitive to the 
new interactions that arise in the presence of these terms. 
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The probable resolution is that the while the full tunneling dynamics is strongly cou- 
pled, the Euclidean calculation gets the correct answer because it correctly captures the 
underlying gravitational thermodynamics of the problem. Thermodynamics tells us that 
the rate is independent of the path chosen, and thus is independent of the fact that Cs may 
depart from unity along the tunneling path. Thus the Euclidean calculation innocently 
gets the correct answer without really capturing the full non-equilibrium dynamics. Since 
gravitational thermodynamics is likely to be fundamental, it seems reasonable to suppose 
that if the contribution from all the higher order noise correlation functions could some- 
how be resummed, the result for the irreversibility ratio TZ would ultimately be consistent 
with the fluctuation theorem, and hence agree with the standard Hawking- Moss result, as 
long as the barrier remains in the regime where the fluctuation determinant can be well 
approximated by a one-loop gaussian integral. 

3.5 Consequences for other instanton tunneling rates 

Another recent calculation of the consequences for tunneling - in the Coleman-de Luccia 
instanton, described in Ref. Q - found a sound speed enhancement to that quantum tun- 
neling rate. As emphasized by [0, though, the Hawking-Moss instanton is the extreme 
case in which tunneling takes place via a maximally broad "bounce", or in a very thick 
tunneling "bubble". The instanton calculation in this case is almost entirely determined 
by the evaluation of the action exactly at the unstable maximum on the tunneling poten- 
tial, where the non-trivial sound speed effects are suppressed. In the opposite extreme, 
Coleman-de Luccia, the quantum tunneling is rapid and proceeds for typically subhorizon 
sized bubbles. Ref. []T9| demonstrates that there are typically a large number of interpo- 
lating instantons between Hawking-Moss and Coleman-de Luccia in the case of canonical 
scalar fields. It is natural to expect that each of these interpolating instantons will also 
exist for DBI and other non-minimal kinetic theories. Though we have not yet pursued the 
question, we expect DBI effects will begin to show up in these other instantons, but will be- 
come progressively smaller as the bubbles get broader and we approach the Hawking-Moss 
instanton. 

4. Summary 

The thermal tunneling rate described by the Hawking-Moss instanton is apparently insen- 
sitive to the strong coupled dynamics of theories with non-minimal kinetic terms as long as 
the barrier lies in the regime where the fluctuation determinant can be well approximated 
by a one-loop gaussian integral. However, the actual non-equilibrium dynamics that a 
tunneling trajectory follows is strongly coupled whenever Cg departs sufficiently from unity 
along the tunneling trajectory. The calculable estimate of the stochastic tunneling rate, 
assuming gaussian noise from subhorizon fluctuations, disagrees with the Hawking-Moss 
calculation in the (cs) ^ 1 case. We have argued that the stochastic tunneling rate and 
the thermal Hawking-Moss tunneling rate should be complementary descriptions of the 
same physics, and consequently it is a surprise that they should appear to disagree. The 
probable resolution of this apparent conflict is that the strongly coupled physics must. 
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when resummed, respect the underlying gravitational thermodynamics that govern de Sit- 
ter space and which are manifest in the local picture. We have explicitly demonstrated that 
the non-gaussian noise fluctuations induced by the subhorizon modes become order unity 
effects during the tunneling process, negating a perturbative calculation of the stochastic 
tunneling rate. Since the fluctuation theorem, the fundamental statement of thermody- 
namics, agrees with the Hawking-Moss instanton, we expect this to be the ultimate answer, 
which implies that the full rate of tunneling is blind to the strongly coupled dynamics of 
the tunneling path. We leave it to future work to hopefully provide a fully non-equilibrium 
derivation of this result and demonstrate the consistency between the local and global 
descriptions of tunneling. 
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